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■ ABSTRACT 

6 

^: 

c/2 . While hundreds of planets have been discovered around field stars, only a few are known in 
■ 

star clusters. To explain the lack of short-period giant planets in globular clusters (GC), such 

^ ■ as 47 Tucane and uj Centauri, it has been suggested that their low metallicities may have 

, prevented planet formation. Alternatively, the high rates of close stellar encounters in these 

in- 

\ clusters may have infiuenced the formation and subsequent evolution of planetary systems. 

I> 
O 
(N 



How common are planets in clusters around normal main-sequence stars? Here we consider 
whether this question can be addressed using data from the Kepler mission. The Kepler field 
of view contains 4 low-density (relative to GCs) open clusters where the metallicities are 
>< : abou. .0.. ,0. eve. M,Ke, a.l .ella. e.eouJe. a.e ,n„eK less .e,uen. in tvpieal 
GCs. We provide detailed A^-body models and show that most planets in Kepler-detectable 
orbits are not significantly perturbed by stellar encounters in these open clusters. We focus 
on the most massive cluster, NGC 6791, which has super-solar metallicity, and find that if 
planets formed in this cluster at the same frequency as observed in the field, Kepler could 
detect 1-20 transiting planets depending on the planet-size distribution and the duration of 
data collection. Due to the large distance to NGC 6791 Kepler will have to search relatively 
faint {Kp < 20) stars for the full extended mission to achieve such a yield. 

Key words: methods: n-body simulations - methods: mimerical - open 
clusters and associations: individual (NGC 6791) - planetary systems 
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2 S. Chatterjee et al. 
1 INTRODUCTION 



Stars are thought to alw ays form in clusters or groups (e.g., ICarpenten 120001 : iLada fc Lada 



2003 



Porras et al. 



20031 ). Most of these clusters dissolve rather quickly under Galactic tides 



and via encounters with giant molecular clouds, forming the apparently isolated stars in 
the field. Multiple independent constraints indicate that our now-is olated star, th e Sun, was 
born in a moderately-sized star cluster with = 10^ - 10^ stars (jAdamsll2010l ). However, 
this understanding leads to an apparent puzzle. Planets within a few AU from the host 



stars are observec. 



le.g- 



Butler et al. 



to 



) e common around normal main sequence (MS) stars in the field 



2006 



Batalha et al. 



20121 ). In contrast, few planets have been discovered 



around MS stars in star clusters. The only planets confirmed in star clusters to date are the 
Jovian-mass plane t in a wide circumbinary orbit around a w hite-dwarf and a millisecond 



pulsar in M4 (e.g.. 



Sigurdsson et al.ll2003l : iFregeau et al.ll2006l ) , the giant planet of minimum 



mass 7.6 M J around e Tauri, a giant star in the open cluster Hyades (jSato et al.l 120071 ). 
and most recently the two hot jupiter planets around F-G stars Pr0201 and Pr0211 in 



A lthough a number of transi t ing planet candidates have been 



reported fe.g.. iMocheiska et al.ll2006 



Lovis &: Mayor 



20071 : 



Miller et al. 



2008 



Pepper et al 



20081 : iBruntt et al.l 120031 ) around clusters stars, the planets recently discovered in Praesepe 



are the only planets discovered around MS stars in a star cluster. Do planets not form around 
cluster stars as often as around field stars? If they do form, does some physical process such 
as stellar encounters destroy planets with detectable orbits? Or have searches simply not 
been sensitive enough? 



^ 10^ Mq) and dense clusters like typical Galactic GCs have low metallicities 
(lHarrislll996 ). Radial velocity (RV) surveys of fie ld stars show that high er metallicity stars 



Massive 



are more likely to host short-period giant planets (IFischer fc Valentill2005l ). Thus, it has been 



suggested that the low metallicities of Galactic GCs may inhibit planet formation around 
stars in these clusters. Nevertheless, recent theoretical studies suggest that planets may be 
able to form even at such low metal licities, especially in short-period orbits {Z ^ O.IZq] 



Kara et al. 



2008 



Johnson Sz Li 



201^ . Alternatively, the high interaction rates in the rela- 



tively crowded regions of star clusters may have led to the destruction of detectable planetary 
orbits (or close encounters may have inhibited their formation via, e.g., disk truncation). In 
contrast to the GCs, open clusters (typically ~ 10^-10'^) have much lower central stellar 
densities. Analytical estimates indicate that in these lower mass clusters stellar encounters 
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should not affect planet formation and the ir subsequent survival, especially for plane ts in 
orbits detectable by transit searches (e.g., iKobayashi fc Idall200ll : lAdams et al.ll2006l ). Of 



course the interaction rate depends on the planet's semimajor axis, and a true semimajor- 
axis-dependent encounter frequency can only be obtained via direct A^-body integration of 
a star cluster with stars hosting planets with a large range of initial semimajor axes. 



Finding a planet around a star in a low- mass cluster is challenging 



hese clusters dissolve 



at a young age under the Galactic tidal forces (iBaumgardt &: Makinoll2003l ). Therefore, the 
observation must be made at a relatively early age before the cluster dissolves. High precision 
Doppler or transit observations are not easy for young stars due to high level of star spots 
and stellar activity. In addition, ongoing star formation, gas expulsion, and presence of high- 
mass (M ^10 Mq) stars in young clusters make these observations hard. Hence, an ideal 
star cluster for a transit search for planets should have the following properties: (1) high 
enough metallicity for planets to form; (2) a sufficiently massive cluster so that it does not 
dissolve too quickly; (3) a sufficiently old cluster so that both star and planet formation 
have been completed. 

Searching for planets in star clusters using transit surveys has been suggested to take 
advantage of the high surface density of target stars (jJanesI Il996l ). Moreover, finding the 
frequency of planets around cluster stars provides us with useful constraints on how often 
planets (within a range of semimajor axes depending on the detection limit) form, mi- 
grate, and survive around a coeval aggregate of stars. Transit surve ys have been used to 



search for hot Jupite rs in 47 Tuc, yielding a significant null res ult ( 



Weldrake et al. 



20051 ) . A similar null result was found for u Cen (jWeldrake et al.ll2008l ). A 



Gilliland et al. 



2000; 



didates around normal 



2006 



2006 



Mocheiska et a. 



Montalto et al. 



2005 



2007 



200q : 



200£ 



e.g., 


Burke et al. 


2006; 


Hartman et al. 


2009; 


Burke et al. 


Miller et al. 


2008 


Bramich et al. 


2005; 


Bramich & Horne 


Rosvick & Robb 


2006 


). Nevertheless, these null results 



from ground based surveys do not yet constrain the planet occ urrence rate in open c 



usters 



201 ih 



to be lower than that expected from observations in the field (Ivan Saders fc Gaudil 
Hence, it would be very interesting to identify a promising candidate cluster where planets 
are expected to be detected. The discovery of planets in such a cluster would provide us 
answers to the apparent puzzle created by the presumed cluster origin of all stars, high 
frequency of short period planets around field stars, and the apparent lack of short period 
planets around cluster stars. Alternatively, a null result from that cluster would potentially 
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put stronger constraints on planet formation efficiency in clusters than are currently present. 
NASA's Kepler mission provides an unprecedented opportunity to investigate this problem 
due to its nearly uninterrupted photometric measurements with high cadence and precision 
over many years. 



Four open clusters lie within the Kepler field of view (FO V), namely NGC 6866 , NGC 6811 



NGC 6819, and NGC 6791, in increasing o rder of their ages (IMeibom et a. 



ters 



l ave solar or higher m etallicities (e.g.. 



2012 



Corsaro et al 



5 X 10^ M 



0) 



Platais et al. 



2011 



Hole et al 



20091 



201I[). A114clus- 



Gunes et al. 



20121) . Amo ng these 4 clusters, NGC 6791 is the most massive (Mt 



Platais et al 



tot 



20111 ) and provides the largest number o 



planets transiting normal MS stars. ]N 



Platais et al. 



2011 



60 MqPc" 



Meibom et al. 



2011 



Platais et al 



GC 6791 is old 



Brogaard et al 



8 ± 1 Gyr; 



targets to search for 



Grundahl et al. 



2008 



2OI2I ). has super-solar metallicity 



([Fe/H] = +0.30: iBoesgaard et al.l l2009l ) and its central density is relatively low (p^ 



2OIII ). much lower than in any GC. The high metallicity indicates 



that planet formation in NGC 6791 should not be inhibited due to lack of metals, and the low 
stellar density relative to typical GCs indicates that stellar encounters are not as common 
as in the GCs. 



In this study we use detailed A^-body simulations with planet-harboring stars and stellar 
binaries incorporating stellar evolution, two-body relaxation. Galactic tidal stripping, dy- 
namical interactions between stars, binaries, and planetary systems, and physical collisions. 
Our goal is to study the semimajor axis-dependent effects of stellar encounters on planetary 
orbits in clusters similar to NGC 6791. Moreover, using initial distributions of planet prop- 
erties based on the Kepler planet candidate list we estimate the number and properties of 
planets that Kepler can detect in NGC 6791 if the planet occurrence rate in this cluster is 
the same as that in the field. 



Our numerical methods are summarized in Section [2l We present our model of NGC 6791 
in Section [3l We then explore various properties indicating the effects of stellar encounters 
on the planetary orbits in our cluster models (Section H]) with particular focus on the model 
that best matches the observed properties of NGC 6791. In Section[5]we investigate whether 
planets in NGC 6791 can be detected by Kepler and discuss the expected properties of these 
detectable planets. We conclude and discuss the implications of our results in Section [61 
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2 NUMERICAL METHODS 

In this section we describe the code we use, how we assign cluster properties (e.g., mass, 
concentration, stellar binary fraction, and fraction of planet-harboring stars; Section [2.ip . 
planetary masses and their orbital properties (Section 12. 2p . We describe how we estimate 
the expected number of Kepler-detectable planets in Section 12.31 



2.1 Modeling of NGC 6791 

We model star clusters with planet hosting stars using a Henon-type Monte Carl o code CMC 



(for c 



2000 



uster Monte Carlo) developed and rigorously tested over the past decade ( 



2001 



Fregeau et al. 



20031: 



Fregeau fc Rasio 



2007 



Chatterjee et al.l 



2010 



Joshi et al 



Umbreit et al. 



2OI2I ). Using this code we can efficiently model realistic clusters with a large number of 
star-systems (A^) and high binary fraction (/;,) including two-body relaxation, binary and 
single stellar evolution, strong encounters (including physical collisions and binary mediated 
scattering interactions), and Galactic tidal stripping at a relatively low computational cost. 
These models can be directly compared with observed star clusters. 

The collective effect of all of the above physical processes over the age of the cluster deter- 
mines the properties of an observed cluster. Although the qualitative effects of each physical 
process is known, the quantitative extent of these effects are hard to estimate without doing 
full simulations. Hence, the initial conditions necessary to create an observed cluster cannot 
be easily inverted. We perform simulations for initial conditions spanning a large grid of 
multi-dimensional parameter space. We focus on modeling NGC 6791 and summarize a col- 
lection of properties for all other models as representative of rich open clusters. We restrict 



the huge parameter space to make calculations tra ctable. For examp 



Galactocentric distance between 5 and lOkpc (e.g.. 



Platais et al. 



e, NGC 6791 has a 



201 ll ). It is also expected 



that the cluster loses about 50 - 70% of its initial mass via Galactic tidal stripping, dynam- 



ical ej ections, and stellar evo . 



[e.g., 



Baumgardt fc Makino 



ution mass loss through win ds and compact object formation 



2003 



Chatterjee et al. 



20101). Considering these a priori con- 



straints we compute a large grid of ~ 200 simulations varying the initial parameters including 
the cluster mass (Md), King concentration parameter {Wq), Galactocentric distance (tqc), 
and virial radius (r^). We compare the simulated cluster properties between 7-9 Gyr with 
the observed properties including the surface number density profile of NGC 6791 to find 
acceptable matches. 
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6 S. Chatterjee et al. 

All our simulated clusters have an initial between 10^ - 10^ and an initial r„ be- 



tween 3-8 pc. The position s and velocities are drawn according to a King model (jKing]|l966 : 



Binney fc Tremaine 



20081 ) with Wq between 3-6. The initial stellar binary fraction (/fe) 



is chosen to be between 0.1 and 0.5. We draw the masses of the stars (or primary stars in 
binaries) from the stellar mass function (MF) presented in KroupJ ( 2001 , Equations 1 and 
2) in the stellar mass range 0.1 — 100 Mq. The mass of each secondary in a stellar binary 
is drawn from a uniform distribution of mass ratios in the range 0.1 - the mass of the 
primary. The initial eccentricities (cj) are chosen to be thermal ( IHeggie &: Hutll2003l ). The 
initial semimajor axis (oj) distribution is flat in logarithmic intervals between 5 times phys- 
ical contact to the local hard-soft boundary of the cluster. The local hard-soft boundary is 
a measure to determine the widest orbit of a stellar binary that will not be disrupted via 
stellar encounters in its cluster environment and depends on the velocity dispersion of the 
stars in the region of the cluster and the binding energy of the binary. The "hard" binaries 
are sufficiently bound to each other and statistically becomes more bound via super-elastic 
stellar encounters. On the other hand "soft" binaries are those that are not sufficiently bound 
and becomes less bound via stellar encounters and eventually gets disrupted via one (or mul- 
tiple) stellar encounters ( iHeggie fc HutI 120031 ). Although all binaries are initially hard, they 
may not remain so throughout the evolution of the cluster. Due to two-body relaxation the 
velocity dispersion near the center of the cluster increases as the core contracts. Moreover, 
binaries sink to the center due to mass segregation where the velocity dispersion is higher 
compared to that at the initial position of a binary. We include this effect in our simula- 
tions and these soft binaries are allowed in the cluster until they disrupt via binary-single or 
binary-binary interactions. The orbital phase angles and orientations are chosen uniformly 
in the full range. 



2.2 Planet properties 

In addition to stellar binaries we also include stars with planetary companions determined 
by a planet fraction /p defined as the ratio of the number of planet host stars (Ap) to the 
number of all star systems (singles or binaries; A^) in the cluster. Each primary can have 
only one companion in our simulations, either a stellar mass companion (we call those stellar 
binaries) or a planetary companion (we call those planetary binaries). Strong interactions 
involving both stellar binaries and planetary binaries in the evolving cluster potential are 
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followed. In any simulation the total fraction of "binaries" (stellar and planetary) is + fp. 
We do not include multiple planet systems or circumbinary (i.e., orbiting a binary star) 
planets due to a limitation in the current version of the code that does not allow us to treat 
hierarchical systems above a binary. 
Planet masses (Mp) are assigned according to a power- law df/dhgMp = M'^-^^ fiHoward et al, 



2OIOI ) between Mp = 1 M0 - 5 Mj where, Mj is the mass of Jupiter. The planetary radii 



(Rp) are assigned a mass-dependent valu e according to Rp = the smaller of -R^ (;^) 



and Rj, where -Rj is the radius of Jupiter (ILissauer et al.ll201ll ). Note that although planets 



with Rp > Rj are observed, we employ the upper limit in Rp to remain conservative in our 
estimated number of Kepler-detectable planets (described later in Section [2^3]) . 

We assign the planetary orbits in two different ways, henceforth Setl and Set2. In Setl 
our focus is to find the effects of stellar encounters on planetary orbits as a function of the 
planetary semimajor axes (op) in a cluster environment. The encounter rate of a binary in a 
star cluster is directly dependent on its semimajor axis a, and is proportional to a^. Hence, to 
estimate the a^-dependent effect of stellar encounters on planetary orbits we need to sample 
Op over a large range. For Setl the initial Op-distribution is flat in logarithmic intervals 
between IQ-^ and 10^ AU. For Setl we use a fixed value of fp = 0.33, a somewhat arbitrary 
choice, but not too far from the overall fraction of Kepler planet candidates observed in the 
field. 

Set2 consists of models with initial conditions resulting in a close match with NGC 6791 
from the large grid of simulations in Setl. Here we focus on estimating the expected number 
of planets detectable by Kepler in NGC 6791 assuming that the initial planet frequency in the 
cluster is the same as observed in the field. In Set2 the orbital period distribut ion for planets 



is gui ded by the observed period distribution of the Kepler planet candidates ( iBatalha et al. 



2OI2I ). Each Kepler planet candidate is weighted by a/R^,, where i?* is the radius of the host 
star, to account for the geometric transit probability. A lognormal distribution for orbital 
period is obtained by fitting the weighted Kepler planet candidate period distribution. We 
draw initial planet orbital periods from this lognormal distribution between 3 and 365 days 
for Set2. We use fp = 45% for Set2. This value of fp is obtained requiring that the rate of 
transiting planets with Rp > 2.5 R^ and period between 3 and 120 days is equal to the rate 
of Kepler planet candidates observed in the field. The Kepler planet candidate list should 
suffer fro m only minor incomp leteness for planets in this range around dwarf host stars with 



Kp < 14 f lBatalha et al.ll2012[ ) 
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In all cases the initial planetary orbits are circular. Since we do not include multiple 



plane ts, there is no excitation of orbits via planet-planet scattering (e.g., I Chatterjee et al. 



20081 ) in our simulations. Thus, any final non-zero e is a result of stellar encounters in the 



cluster. In (Section l4.ip we argue that indirect planet-planet instabilities, i.e., planet-planet 
scattering triggered by stellar fiy-bys, should be a small percentage effect in clusters like 
NGC 6791. The orbital phase and orientation angles are drawn uniformly in the full range 
in both sets. Planetary companions are assigned only around host stars with ^ 2 Mq in 
both sets since transiting planet searches focus on these systems. 



2.3 Detectability of Planet Transits by Kepler 



r simulations we 


;rack 5 


Chatterjee et al. 


2010) 



(e.g., Idiatterjee et al.l 120101 ). Thus, at any given time in the evolution, we can extract the 
stellar radii (i?^,) and bolometric luminosities (L^) of planet host stars and the orbital prop- 
erties of the planets. The bolometric luminosities are converted to standard colors using the 
standard filter band passes (Johnson i?, V) and synthetic stellar spectra obtain ed from stellar 



atmos pheric models dependent on the stellar metallicities, and surface gravity (ILejeune et al. 



1997a 



b|). The -B, V band magnitud es are then converte d to SDSS g, and r magnitudes using 



the transformation equations from 



Smith et al. 



(I2OO2I ). The Kepler magnitude Kp is then 



calculated using the standard Kepler conversions obtained from the Guest Observer (GO) 



website of Kepler (http://keplergo.arc.nasa.gov/CalibrationZeropoint.shtml). Note that the 
conversion e q uatio ns in the Kepler GO website is slightly different from those given in 



Brown et al. 



(I2OIII ). We find that the Kp values calculated using the conversion equations 
in the Kepler GO website match better with the actual Kp values in the Kepler target list 
(lower rms d i fferen ce) compared to the Kp values calculated using the equations given in 
Brown et all J201l[ ). 

Whether or not a planet would be detected by Kepler depends on the ratio of the transit 
signal strength and the noise metric from all contributing factors (including photon noise 
and intrinsic stellar variability) for Kepler in a specified temporal length for the host. For a 
noise metric we use the combined differential photometric precision (CDPP) measured over 
6.5 hours. We assign the CDPP values in ppm for the planet host stars in our simulations 
using a polynomial fit: 

log CDPP{Kp) = 10.488 - 2.mSKp + 0.341827^^^ - 0.0170457^^^ + 0.000331467^/ (1) 
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to the CDPP values measured by Kepler and given in 



Gilliland et al 



( 1201 ll . their Figure 4). 



Note that in general extrapolating the CDPP values to large Kp values (Kp = 20 n eeded in 



Gilliland et al. 



201 If ) can be 



our study) using CDPP data in the range Kp < 16 (given in 
dangerous. However, an investigation by Gilliland (private communication with the referee 
Gilliland) unveils that the CDPP values predicted using this extrapolation equation indeed 
provides CDPP values in acceptable agreement with real data even for Kp = 20. The values 
predicted by the extrapolation equation provides a more conservative estimate compared to 
the real data (private communication with the referee, Gilliland). The signal to noise ratio 
(SNR) of a transiting planet is calculated for each planet using 

rnpp U .^hr/ 



SNR 



(IHoward et al. 



(2) 



2OI2I ). where is the number of transit events within a given length of data 



collection by Kepler (e.g., 1, 3.5, and Syr are considered), and tdur is the duration of transit, 
given by 

a 



t. 



dur 



2i?. 



GAL 



0.5 



(3) 



assuming a circular orbit and central transit. If the SNR > 7, we consider the planetary orbit 
to be "detectable" by Kepler. Each detectable planet is then weighted by the geometrical 
transit probability Ptr = -R^/o to estimate the actual number of transiting planets detected 
by a single observer. 

The predicted number of detectable transiting planets can vary due to statistical fluc- 
tuations in their periods, sizes, and which planets were assigned to which stars in our sim- 
ulations. Therefore, we vary the seed of the random number generator and generate 4 re- 
alizations of the initial conditions for each set of cluster input parameters that gives us a 
good match for NGC 6791 to estimate the size of the statistical fluct uations. Note that th e 



estimated CDPP values for a given Kp is obtained for stars on the MS ( 1 Gilliland et al. 



20111). 



If the CDPP values are significantly different around giant stars, then the noise metric esti- 
mated using Equation [1] may be inadequate. Nevertheless, we will show that the majority of 
Kepler detectable transiting planets are around bright MS stars (and some subgiant stars), 
hence, this caveat is not expected to be serious for the purpose of this study. 



3 SIMULATED MODEL OF NGC 6791 
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Table 1: List of properties for our simulated model of 
NGC 6791. The initial parameters are given as well as 
the properties of the model cluster at 8 Gyr. Two columns 
for the final properties are based on including all stars in 
the model, and a subset of stars sa tisfying q < 22. The 
observed properties are taken from 
for comparison. 



Platais et al. 



(120 111 ) 



Cluster Property 


Initial Value 


Final Value 


Observed Value 




All Star g <22 


g<22 


N 


8 X 10^ 


1.8 X 10^ 3.8 X 10^ 




Total Cluster Mass (M©) 


6 X 10^ 


10^ 4.5 X 10^ 


5 X 10^ 


Concentration parameter (log ^) 


1.1 


0.91 


0.9 


Virial Radius 


8 pc 






Galactocentric Distance (kpc) 


8.5 


8.5 


5-10 


Stellar binary fraction (ft,) 


0.30 


0.35 




Fraction of planet-harboring stars (/p) 


0.33 


0.27 





To compare our models from Setl with the observed NGC 6791 we need to "observe" 
our cluster model. We compare the surface number density profiles from o ur models, a 
basic observable for clusters, with that of NGC 6791 at tci = 8 ± 1 Gyr ( jPlatais et al. 



20111 ). In order to compare with the observed cluster we convert the bolometric luminosities 



obtained directly from CMC to standard colors using the standard filter band passes and 
synthetic stellar spectra obtained from stellar atmosph eric models dependent on the stellar 



metallicities, and surface gravity ( iLejeune et al. 



1997aU bl). We only include objects brighter 



than < 22 to create the surface density profiles. A good match in the surface brightness 
profile indicates that the total luminosity (a proxy for mass) of the cluster as well as its 
radial distribution is well modeled. Indeed, in observed clusters the surface density profile is 
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Figure 1. Surface density profile for the observed cluster (black dots) and the best- match simulated cluster at three times 
between = 8 ± 1 Gyr (blue, green, and brown in increasing order of td) during its evolution. For the simulated cluster only 
stars with g < 22 are considered to be compa tible with the obse r ved d ata. The observed data is obtained using the Dexter 
data extraction applet l lDemleitner et al]|200ll ) from lPlatais et al] l l201ll ). One of the key questions for this study is how much 
planetary orbits are perturbed in a cluster environment. This is directly dependent on the density profile of the cluster. The 
close similarity in the density profiles for the observed NGC 6791 and our theoretical models indicate that the dynamical 
interaction rate in our model should be very similar to that in NGC 6791. 



often used to estimate many global structural properties including central density (pc), core 
radius (rc), and half-light radius (rhi). 

Note that for this study we do not compare various stellar populations. Clusters are 
birth places for a number of exotic stellar populations including blue straggler stars (which 
should be created mostly via mass-transfer in a stellar binary giv en the low central density of 
NGC 6791), and low-mass X-ray binaries (jPooley fc Hutll2006l ). These populations depend 
closely on the details of the initial stellar MF, stellar binary properties, the stellar density 
profile, and the dynamical evolution of the cluster. While numerical modeling of star clusters 
focusing on reproduction of these individual exotic stellar populations is a very interesting 
and active area of research, it is beyond the scope of this study. For our purposes, the most 
important aspect in a cluster environment is its density distribution since the stellar density 
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Figure 2. A synthetic CMD for the best-match cluster model for NGC 6791. The single and binary MSs are clearly seen as 
well as the red giant branch. We find that the best chance for detecting planets is around stars residing near the high end of 
the MS (box). The giants are brighter, but they are also much bigger (the radius can be 10'^ times compare d to the radius of 
the star during MS) resulting in a drop in the signal to noise ratio for Kepler. The number of blue stragglers ( Sandagc 1953) in 
our model is lower compared to what is observed in NGC 6791. The blue stragglers in NGC 6791 are expected to be created via 
mass transfer in binaries due to the low central density in the cluster and should not be indicative of the dynamical interaction 
rate which is dependent on the density profile of the cluster (Figure [TJ. 



distribution directly affects the interaction rates and in turn affects how the properties of 
planetary orbits are expected to change in the cluster. 

The initial conditions that produce a good match with NGC 6791 are presented in Table 
[U Table [T] also lists the final properties of the model cluster at age t = 8 Gyr and some 
observed properties of NGC 6791. Figure [1] shows the comparison between the observed 
surface density profile of NGC 6791 and that for our simulated best-match model at three 
different snapshots between tci = 8 ± 1 Gyr. 

By tracking both the dynamical evolution of the cluster objects in the overall cluster 
potential and the stellar and binary properties as they evolve we can create realistic synthetic 
color-magnitude diagrams (CMD) for our simulated models. Figure |2] shows an example of 
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synthetic CMD created for the cluster model showing best match for the observed NGC 6791 
at tci = 8 Gyr. The single and binary MS s, the MS turn-off, and the Giant branch are clearly 

d (ab out 8 Gyr; 



visible. NGC 6791 is o 

us (4kpc; 



Platais et al. 



Grundahl et al 



2008[ ) and it is relatively far from 



201l[ ). Hence, the MS stars in this cluster are not very bright. Note 



that the single MS turn-off is at ~ 17 Kp and the binary MS turn-off is near 16i^p. 

We summarize some key properties of all our simulated clusters in the wide multidimen- 
sional grid from our search in Table [2] available in its entirety in the electronic version. The 
best-match model for NGC 6791 is listed as runl67 in the Table [2J In the rest of this study 
we will focus on the results using these cluster parameters. 



© 2002 RAS, MNRAS OOO.nHiTI 



14 S. Chatterjee et al. 

Table 2: List of properties for all our simulated cluster 
models at different cluster ages {td). This table is pre- 
sented in its entirety in the online journal. A portion 
is presented here for guidance to its form and content. 
The initial parameters for each model including the King 
concentration parameter Wq, virial radius r^,, Galactocen- 
tric distance tgc are given below each run name. Vari- 
ous properties are listed at a few t^i values. The total 
cluster mass M, total cluster object number N , central 
density pc? stellar binary fraction /;,, fraction of planet 
host stars fp are listed. The ratio of the number of free- 
floating planets bound to the cluster {Npjf abound) to the 
number of planets bound to their host stars in long or- 
bits Np{a > 10 AU) is denoted as fpjj^bound- The ratio 
between Np and the number of planets escaped from the 
cluster potential and still bound to their hosts Np^^sc is de- 
noted by fp^esc- The ratio between Np^^sc and the number 
of free-floating planets escaped from the cluster potential 
{Npjf^esc) is denoted by fpjj^esc- Cluster models that dis- 
solve due to Galactic tides within a Gyr because of their 
low mass and compactness are excluded from the list. 



Name 


td(Gyr) 


M (10^ Mo) 


A^/10=^ 


Pc (M0PC-3) 


fb 


fp 


fp,ff,bound 


fp,esc 


fp,ff,esc 


runl 


0.00 


13.4 


20.0 


392.7 


0.10 


0.33 


0.0006 







H^o = 3 


0.27 


9.1 


19.2 


13.7 


0.10 


0.33 


0.001 


0.03 





= 5.0 pc 


0.55 


8.3 


18.0 


11.4 


0.10 


0.33 


0.002 


0.1 


0.002 


Tgc = 7 kpc 


0.82 


7.6 


16.8 


10.7 


0.10 


0.33 


0.002 


0.2 


0.006 




1.09 


6.9 


15.2 


10.1 


0.10 


0.33 


0.004 


0.3 


0.005 




1.36 


6.3 


13.6 


9.8 


0.10 


0.33 


0.004 


0.5 


0.005 




1.64 


5.5 


11.6 


7.4 


0.11 


0.32 


0.002 


0.7 


0.005 




1.90 


4.7 


9.6 


8.0 


0.11 


0.32 


0.005 


1 


0.005 



Continued on Next Page. . . 
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Name 


tci Gyr 


M (10^ Mo) 


A^/10=^ 


Pc (Mepc-^) 


fb 


fp 




f 


IT) f f P 'iC 




2.17 


3.8 


7.3 


9.0 


0.12 


0.32 


0.005 


2 


0.004 




2.45 


2.9 


5.3 


9.0 


0.12 


0.31 


0.009 


3 


0.004 




Z. iZ 


Z.Z 


o.o 


y.o 


U.io 


U.di 


U.Ui 





U.UU4 


run2 


0.00 


27.1 


40.0 


657.9 


0.10 


0.33 










Wo = 3 


0.34 


18.4 


39.3 


27.4 


0.10 


0.33 


0.005 


0.01 


0.07 


= 5.0 pc 


0.68 


17.4 


38.6 


21.0 


0.10 


0.33 


0.009 


0.03 


0.05 


roc = 7 kpc 


1.03 


16.6 


37.6 


18.7 


0.10 


0.33 


0.01 


0.06 


0.03 



4 CLUSTER EFFECTS ON PLANETS 

Here we investigate the effects of stellar encounters on planetary orbits as a function of the 
planetary semimajor axis in open clusters similar to NGC 6791. 

4.1 Stellar encounter and planetary orbits 

Figure [3] shows the ap-distributions of the planetary objects in our numerical model of 
NGC 6791 (Table [1]) at t = and 8 Gyr. The Op-distribution at 8 Gyr remains very similar 
to the initial distribution apart from a moderate spreading of the Op ranges. For Op ^ 0.02 AU 
the planetary orbits can shrink due to tidal damping. Moreover, stellar encounters make hard 
planetary orbits harder. Similarly, soft orbits expand via stellar encounters. In addition, 
stellar-evolution driven mass loss from the host star also expands planetary orbits. However, 
the overall shape of the ap-distribution remains more or less unchanged indicating that stellar 
encounters in this cluster do not significantly change most of the planetary semimajor axes 
with Op ^ 100 AU. 

A more direct way to probe the importance or the lack thereof of stellar encounters in 
changing the planetary orbits in clusters similar to NGC 6791 is investigating the number of 
planets with non-zero eccentricities as a function of Op. We initialize all planets on circular 
orbits and all non-zero eccentricities are results of stellar encounters. Figure H] shows the 
ap-distribution of planet orbits with e ^ some given value. Overall, only about 10% of all 
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Figure 3. Initial (black, solid) and final (t = 8 Gyr; red, dashed) probability density function of the planetary semimajor axis 
for all planets (ap) in the cluster. Note that the fraction of planetary orbits with very small ^ 0.02 AU gets reduced. This is 
mainly due to some stars evolving and engulfing close-in planets, and shrinking of orbits due to tides. A small fraction of large 
ap-orbits expand via interactions as well as stellar evolution driven mass loss. Nevertheless, the shape of the distribution for 
most orbits is unchanged indicating low rate of interactions. 



planetary orbits acquire e > 0.1 via interactions. Furthermore, the close-in orbits (op < 
10 AU) are mostly unperturbed. Since transit surveys preferentially find planets with small 
semimajor axes, transit-detectable orbits are expected to remain unchanged against possible 
stellar encounters in clusters like NGC 6791. Thus, these orbits are not expected to be 
different from those in the field. We find the same general result for all our simulated models 
in Setl although the exact fraction of significantly perturbed planets depends on the stellar 
density of the clusters (see Table [2]). 

Our simulations include only single planet systems. In multi-planet systems indirect in- 
stabilities may occur when a stellar encounter changes the outer pla net's orbit and this ex- 
citation in turn excites t he orbits of the inner planetary system (e.g., 



Zakamska fc Tremaine 



2004 



Boley et al. 



20121 ). Nevertheless, the low encounter rates even for large-Op orbits 

© 2002 RAS, MNRAS OOO.ITMTI 



Planets in clusters 17 



10" 



lO' 



10" 



e~0 (0.90) 
e>=0.2 (0.04) 
e>=0.4 (0.03) 



-1 
log (a^/AU) 



— e>=0.6 (0.03) 

— e>=0.8(0.02) 





Figure 4. Semimajor axis distributions of all planets with e ^ the given value. Since all planet orbits are initially circular a non- 
zero e denotes that the planet have interacted via stellar fly- by, or more complicated binary-binary interactions. Close orbits, 
detectable by transit searches remain mostly unperturbed in a cluster like NGC 6791. Only planets with large-Op (ap ^ 10 AU) 
orbits have a significant chance of being perturbed. The fraction of planetary orbits contained inside each histogram is listed 
in parentheses. Only a small (~ 10%) fraction of planetary orbits gain e via encounters. 



(op ^ 10 AU) indicate that such indirect instabihties should be rare and hmited to a few 
percent effect in clusters similar to NGC 6791 (Figure S]). Of course, independent of stellar 
enco unters, planet-planet in t eractions in mult i-planet systems can perturb planetary orbits 



(e.g., 



Chatterjee et al 



2008 



Naozetal 



201l[ ). Indeed, planet-planet interactions is a com- 



mon source of perturbations for planetary systems orbiting field stars and cluster stars and 
the excitations due to stellar encounters serve as an additional source of perturbation in 
star clusters. Our simulations suggest that for clusters similar to NGC 6791, perturbations 
from stellar encounters should be a sufficiently small effect so that planets detectable via 
the transit method should not be much different from those in the field. 

Even the very long-period planets (a ~ 10 - 100 AU) mostly remain bound to the host 
stars for the whole lifetime of clusters similar to NGC 6791 (FigureE]). The largest orbits have 
the highest interaction cross-section. These orbits are dynamically soft, resulting in further 
expansion of these orbits via typical stellar encounters in these clusters. In addition, stellar 
evolution driven mass loss further expands some of these orbits. Our simulations suggest 

© 2002 RAS, MNRAS OOO.nHiTI 



18 S. Chatterjee et al. 




Figure 5. Histogram for the fraction of free-floating planets (/p,//, normalized so that the integrated area under the curve is 
unity) for our collection of simulated cluster models /p,// = {Np -f-f i,aund + ^p,f f,e3c)/^p{o,p > 10AU,t = 0), where Np -f-f is 
the number of free-floating planets still bound to the cluster potential, Np jj ,,^^ is the number of free-floating planets escaped 
from the cluster potential via ejections and tidal truncation due to Galactic tides, and Np(ap > 10AU,t = 0) is the initial 
number of planets with wide (a > 10 AU) orbits. The solid (black) and dashed (red) lines are for all models and models similar 
in properties with those of the NGC 6791, respectively. The data from all simulations are listed in Table [2] 



that a small fraction of planetary orbits can expand to Op ~ 1000 AU and still remain bound 
to their hosts. When such clusters dissolve near the Sun, planets with these large-a orbits 
would populate the field and provide potential targets for direct imaging planet searches 
while the planets are still young and luminous. 



4.2 Free-floating planets 

Only a small fraction of planetary orbits are significantly perturbed (Figure [3]). Most of 
these planets initially have large-a^ (> 10 AU) orbits (Figure H]). A fraction of these large-Op 
orbits can be disrupted via stellar encounters creating free-floating planets. Mass segregation 
drives the free-floating planets away from the cluster center, but the free-floating planets can 
remain bound to the cluster p otential for many relaxation times near the outskirts of the 
cluster ( iHurley fc Sharal l2002[ ) . In our simulations we find that these free-floating planets 
reside near the tidal boundary of the cluster confirming previous results. Some free-floating 
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planetary objects 




t (Myr) 

Figure 6. Cumulative fraction of planetary objects leaving the cluster potential as a function of cluster age. Planets may 
leave the cluster due to dynamical ejections resulting from strong binary-mediated interactions predominantly in the core of 
the cluster or due to Galactic tidal stripping from the tidal boundary of the cluster as a function of cluster age. Solid (red), 
dashed (black), and dotted (blue) lines denote all planetary objects, planetary objects still bound to a host star, and host less 
planetary objects, respectively. All numbers are normalized by the initial number of planets N{0). Note that only a few planets 
lose their host stars due to stellar interactions in NGC 6791. Most planets are lost from the cluster as the planets' hosts escape 
from the cluster potential. 



planets can remain bound to the cluster for a few billion years after formation. Since the 
free-floating planets leave the core of the cluster due to mass segregation quickly, further 
interaction of free-floating planets with other planetary orbits is unlikely. 

The number of free-floating planets generated by stellar encounters remains low 100), 
a direct effect of low encounter rates in these clusters. This number depends on the actual 
number of planets with sufficiently wide orbits for interaction with other stars or binaries 
in the cluster. Hence we define a quantity fpjf, bound = Npj j abound /Np{a > 10 AU), where A^j, 
is the number of planets bound to hosts in the cluster, Np{a > 10 AU) the same but with 
orbits wider than 10 AU, and Npjf abound is the number of free-floating planets bound to the 
cluster potential. For clusters similar to NGC 6791 fpjf, bound remains a few percent for all 
times. In denser clusters in our simulations fpjj, bound can be as high as ~ 20%. At td = 
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fpjf, bound = by construction. In general, fpjf abound grows until an age td ~ 1 Gyr followed 
by a steady decrease as the cluster grows older due to preferential stripping of free-floating 
planets via Galactic tides (Table 

Microlensing surve ys indicate that th ere is a large number of free-floating (or wide orbit) 



planets in the Galaxy (jSumi et al. 



201ll ). Microlensing surveys have also detected planetary 



mass candidates in the young open clusters a Orionis and Orion (IZapatero Osorio et al. 



2000. 


2002; 


Bihain et al. 


2009) 



20091 ). Hence, it is interesting to investigate what fraction of wide 



planetary orbits are disrupted to create free-floating planets that may or may not remain 
bound to the cluster. For this purpose we define a quantity fpjf = Npjf/Np^i{ap > 10 AU), 
where Npjj is the total number of free-floating planets (bound to and escaped from the 
cluster), and Np^i{ap > 10 AU) is the initial number of wide planetary orbits. Figure |5] shows 
the distribution for /p,// for all our simulated cluster models and for all cluster ages. The 
values of fpjf are typically about a few percent for our simulated models. Thus, if stellar 
encounters dominate production of free-floating planets then most microlensing planets are 
likely wide orbit planets, rather than free-floating planets. However, note that this is a lower 
limit for /p,// for clusters similar to our simulated models. The actual number of free-floating 
planets could be significantly higher due to the cori tribution from planet-planet scattering 
in multi-planet systems (e.g.. IChatterjee et al.ll2008l ). 



4.3 Number of planets and cluster age 

The total number of star systems (single and stellar or planetary binaries) bound to the 
cluster steadily decreases over time due to Galactic tidal stripping. As a result the total 
number of planets in the cluster also steadily decreases over time. However, the ratio of 
planet host stars in the cluster to all star systems in the cluster {fp = Np/N) remains 
roughly constant (Table [2]). For dissolving clusters the ratio fp^esc = Np^esc/Np grows with 
time, where, Np^esc is the number of planets escaped from the cluster potential. In our 
simulations, /p,esc can become very large depending on the cluster age and properties (~ 100 
e.g., runl45). For a completely dissolved cluster the limiting value is fp^esc oo. For our 
best-match model of NGC 6791 fp^esc = 4 at a cluster age td = 8 Gyr (runl67. Table El). 
Thus, for every planet in this cluster there are 4 planets that are now in the field due to the 
slow dissolution of the cluster. 

This decrease in Np is dominated by the host being lost from the cluster. The number 
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of planets lost from the cluster still bound to their host stars is typically ^ 10^ times the 
number of free-floating planets lost from the cluster at all cluster ages (Figure IH])- This is 
simply because complete ionization of planetary orbits is not common in clusters such as 
NGC 6791 even for the broad range in ap considered in our simulations in Setl. The fraction 
of free-floating escaped planets with respect to escaped planets still bound to their hosts 
if p J f, esc) are given in Table [2] for all our simulated clusters at various ages. Note again that 
the fraction of free-floating planets in the field presented here is a lower limit based solely on 
the stellar interactions in the cluster of origin for the planet-host stars and does not account 
for possible contribution from planet-planet scattering. 



5 DETECTABILITY OF PLANETS IN NGC 6791 USING KEPLER 



We have established that stellar dynamics in clusters similar to NGC 6791 alone cannot 
significantly alter small-Op planetary orbits that are detectable via transit surveys like Kepler. 
In addition, we kno w that NGC 6791, already in the Kepler FOV, has super-solar metallicity 



Platais et al 



201ll ). Hence, formation of planets should not be reduced due to lack of solids 



in protoplanetary disks in this cluster. Now we use the simulations from Set2 (Section [2]) 
to estimate how many planets Kepler could discover in NGC 6791 assuming that planets 
form around cluster stars at the same frequency as observed in the field by Kepler and that 
their initial planetary orbits have similar properties to those orbiting around field stars. 
Note that initially we do not in clude any boost in plane t occurrence rate due to the super- 



solar metallicity of NGC 6791 ( IFischer fc Valentil l2005l ). We discuss briefly the estimated 



effects of the high metallicity of NGC 6791 on the predicted number of planet detections 
later in Section 15.21 The initial conditions for our models in Set2 are taken from the best- 
match model in Setl (Tabled]). Set2 consists of 4 different realizations of the same initial 
cluster. The period distribution of the planets follow the di stribution seen in t he Kepler 
planet candidate list after de-biasing for transit probability f Batalha et al. 2OI2I ). Initially 
we consider a planet mass distribution estimated from the RV surveys. However, later in 
this section (Section 15. 2p we explore how the results would change if we adopt other planet 
size distributions including ones based on Kepler observations. The fp is chosen based on 
the observed field planet occurrence rate (Section l2.2p . 

Figure [7] shows a cumulative histogram of the number of transiting planets {Np{< Kp)) 
detectable by Kepler as a function of the Kepler magnitude. Kp = 16.5 is near the MS turn- 
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V 




Figure 7. Cumulative histogram of the number of Kepler-detectable (SNR > 7) transiting planets, A^j 
of Kp. Solid (black), dashed (red), and dotted (blue) lines assume Kepler missions extending for 1, 3.5, 
The black, red, and blue shaded regions denote the range of Np(< Kp) between all 4 realizations of 
same initial conditions indicating statistical fluctuations. Note that Np{< Kp) rises sharply near the MS 
at Kp = 16.5 (Figure [2]l. About 10 planets are expected to be detectable by Kepler by analyzing just 
relatively fainter stars (16 < Kp < 20) are also fully analyzed. If an extended mission provides 8 years of 
of detected planets could grow to about 20 including several smaller planets. Note that these numbers 
overall normalization fp and also the adopted planet-size distribution (see text; Table 



, (< Kp), as a function 
and Syr, respectively, 
simulations using the 
turn-off for NGC 6791 
1 yr of Kepler data if 
data, then the number 
are dependent on the 



off of NGC 6791 (Figure[2]) and the steep rise in A'p(< Kp) near Kp = 16.5 is reflective of the 
stellar properties in the cluster. Although the giants {Kp < 16.5) are much brighter than the 
MS stars, they are less numerous. In addition, the giants also have a much larger R^. Hence 
the transit SNR is not sufficient for most planets to be detected around these giant stars. 
However, there is a chance that Kepler may detect a few giant planets around low-luminosity 
giant stars. Most planet detections in NGC 6791 with Kepler will be around high-luminosity 
MS star hosts (20 > Kp > 16.5). Our results suggest that if planet occurrence rate in 
NGC 6791 is similar to that observed in the field, Kepler could detect a few to 10 planets 
depending on the adopted planet-size distribution even with a single year of data collection. 
A total of about 3-15 and 4-20 planets are expected to be detected from 3.5 yr and 8 yr of 
data collection by Kepler. These ranges reflect uncertainties from the planet-size distribution 
(Table [3]). Relatively faint stars {Kp < 20) must also be analyzed to achieve this yield. The 
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V 




Figure 8. Cumulative histogram of the number of Kepler-detectable (SNR > 7) transiting planets as a function of the 
planetary radius. The lines, colors and shades have the same meaning as in Figure [7] For observations over one year, the 
median of the radii of detected planets is Rp 7 Ris- An extended mission is expected to detect several smaller planets. An 
extended mission with 8 years of data collection may result in a handful of planets with Rp Ri 2 ijgj . Note that the actual 
number will vary depending on fp and the planet-size distribution. 



lack of detections of Kepler candidates in this cluster to date may be due to inadequate 
analysis of these faint target stars. 

Figure |8] shows the cumulative number of transiting planets detectable by Kepler in 
NGC 6791 as a function of the planetary radius (-Rp). About 35% and 85% of all Kepler- 
detectable planets are expected to be larger than Saturn and Neptune, respectively, if 1 year 
of data is analyzed. Increasing the duration of analyzed Kepler data increases the number 
of detections for smaller planets. For example, the median Rp for the predicted number of 
Kepler-detectable planets is about 7, 5, and 4/?® when analyzing 1, 3.5, and Syr of Kepler 
data, respectively. The expected numbers of detectable planets with Rp ^ -RNeptune are about 
2, 5, and 9 when analyzing 1, 3.5, and Syr of Kepler data, respectively. Note that the actual 
numbers of detections may be higher or lower by a few planets due to statistical fluctuations 
as shown by the shaded regions in Figures [7] and [HI The sharp increase in the detectable 
planet numbers in the Rp = Rj bin is reflective of our assumed mass-to-radius relation which 
does not allow our planets to have Rp > Rj (Section |2]). Without this assumption some of 
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log(r/PC) 



Figure 9. Cumulative histogram of the number of the Kepler-detectable planets (SNR > 7) as a function of the planet's 
sky-projected distance from the center of the cluster. Curves and shades have the same meaning as in Figure [7] About 70% 
of all planets reside outside the core {tc = 3.8 pc ^ 3.3') of the cluster. A little over 60% of these planets reside outside 
the half-light radius (r^j = 5.1 pc fa 4.4') of the cluster. The vertical thin dotted lines mark the projected distances from 
the cluster center where the number of stars within one EE95 for Kepler drops below 1 assuming an EE95 size of either 
4.5 or 5.6 pixels (rEE95 = 5.6, 7 pc « 4.8', 6', respectively). Although the actual numbers can vary depending on fp and the 
planet-size distribution, the overall radial distribution should not change. The vertical thick dashed lines denote the boundary 
of the 200 X 200 square pixel block centered on the cluster center for the superaperture where all pixel level data is being 
downloaded. These two lines (for superaperture) denote radial distances from the center of the cluster along a side (6.7') and 
along a diago nal (9.4') of the squ are superaperture box. For reference, one PC is equivalent to 1.16 arcsccond at the distance 
of NGC 6791 l lPlatais et al.ll201ll) . 



the planets in this bin would spread into bins with larger radii. Almost all detectable planets 
for 1, 3.5, and Syr observations are expected to be in orbits with periods below about 40, 
80, and 100 days, respectively. 

Figure [9] shows the cumulative number of Kepler-detectable planets as a function of 
the sky-projected distance from the center of the cluster. Most Kepler-detectable planets 
are not in the most concentrated part of the cluster near the center. In fact, only about 
30% of detectable planets are expected to reside inside the projected Tc = 3.8 pc ~ 3.3' of 
the cluster. A little over 60% of the detectable planets are expected to reside outside the 
projected half-light radius r^i = 5.1 pc ^ 4.4' of the cluster. This is due to mass segregation 
in the cluster and does not include observational selection effects due to crowding. 

In general, crowding of stars can become a potential problem for transit searches in a 
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star cluster. The extent of crowding for each star in the Kepler FOV depends on the number 
of pixels enclosing 95% of the star's flux (EE95). The EE95 values for the module contain- 
ing NGC 6791 in the Kepler FOV varies between 4.5 - 5.6 pixels (obtained from the Ke- 
pler Instrument Handbook; http: / /keplergo. arc. nasa.gov/calibration/KSCI-19033-001.pdf) 
depending on the quarter of observation. Adopting EE95 of 6 pixels and an angular scale per 
pixel of 3.98 arcsec onds (equivalent to an area of about 0.2 pc^ at the distance of NGC 6791; 



Bryson et al. 



20101), the number density of stars with Kp < 22 in NGC 6791 is < 1 at a 



cluster-centric two dimensional {2D) projected distance (rgEgs) of about 7pc 5.6 pc if 
EE95 is 4.5 pixels; Figure H]). Between 33 - 50% of all otherwise detectable planets in our 
simulations reside outside rsEgs = 5.6 - 7 pc ~ 4.8' - 6'. Even at the center of NGC 6791 
the number of stars per pixel is < 1. Hence, crowding should not be a fundamental limit for 
planet search in NGC 6791 using Kepler. However, the current pipeline may not be sufficient 
due to technical issues associated with targeting and analysis of stars in a dense field. 

The above estimated numbers of Kepler-detectable planets in NGC 6791 can change 
based on various parameters including the intrinsic planet frequency, the planet-size distri- 
bution, metallicity of NGC 6791, limited data availability and crowding. The effect of the 
assumed intrinsic planet-occurrence frequency, denoted by fp in this study, on the number of 
detectable transits is straightforward. If all other parameters such as the period distribution 
and the planet-size distribution are kept fixed, the number of transit detections by Kepler 
simply scales linearly with fp. In the following sections we discuss the effects of the other 
parameters on the number of detectable transits by Kepler in NGC 6791. 



5.1 Effects of planet-size distribution 

For a fixed fp, estimating the change in the number of Kepler detectable planets depending 
on the distribution of planet sizes is relatively more complicated. Assuming that the planet 
size distribution is described by a simple power-law of the form df/dRp oc Rp, we estimate 
the dependence of the number of expected transit detections by Kepler in NGC 6791 as a 
function of the power-law exponent a in the following way. We extract the planet-harboring 
stars in our simulations at 8 Gyr but now change the planetary sizes based on several different 
power-law distributions. For each value of a, 200 realizations are generated to estimate the 
statistical fluctuations. Each of these planetary systems are then analyzed using the same 
method described in Section 12.31 We use a range of a values spanning those reported in 

© 2002 RAS, MNRAS OOO.nHiTI 



26 S. Chatterjee et al. 

literature. We report the expected number of transit detections in NGC 6791 by Kepler as 
a function of a in Table El Depending on the value of a, the expected number of detections 
can have values over an order of magnitude. 

Table 3: The expected number of Kepler-detectable tran- 
siting planets as a function of a, where a is the power- 
law exponent for an assumed planetary size distribution 
of the form df/dRp cx R^. Sources (where exist) for a 
values are given in the "Comments" column. These esti- 
mates do not include any boost of planet occurrence due 
to the super-solar metallicity of NGC 6791. These esti- 
mates assume that data from the full cluster is available. 



a 


Number of Detected Transits 


Comments 




1 year 3.5 year 8 year 





-1.5 


12 ±2 


16 ±2 


19 ±2 




-1.99 


11 ±2 


15 ±3 


19 ±4 


these simulations: Howard et al. f2010l 


-2.7 


3 ±2 


6 ±2 


8 ± 2 


Youdin f2011) 


-2.92 


2 ± 1 


4±2 


7±2 


Howard et al. f2012) 


-3.5 


1 ± 1 


3 ±2 


4±2 





5.2 Effects of metallicity 



The p lanet occurrence rate varies with the metallicity of the target stars ([Fischer fc Valenti 



20091). 



20051 ). NGC 6791 has a super-solar metallicity ([Fe/H] = +0.30; iBoesgaard et al. 
Analysis up to this point did not account for the fact that the high metallicity of NGC 6791 
is likely to boost planet occurrence in this cluster relative to the Kepler target list. The 



corre 
(e.g., 



ation between planet occurrence and metallicity is well establish ed for giant planets 



Fischer &: Valenti 



2005 



Valenti Sz Fischer 



for planets smaller than 4 (e.g. 



Sousa et al 



2008 



2011 



Santos et al. 



Buchhave et al 



20111) bu t is not as strong 



20121) . Here we explore 
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how much the super-solar metaUicity of NGC 6791 could change the estimated number of 
detectable transiting planets by Kepler in this cluster. To remain conservative we assume 
that the planet frequency is increased due to metaUicity only for planets larger than Saturn. 
We further assume that the distributions of orbital properties (e.g., period distribution) of 
hot jupiters are independent of t he metaUicity. 



Based on the RV observations 



Fischer fc Valentil (120051 ) fit for the planet occurrence rate 



and the metaUicity of the host stars and find 
U{[Fe/H]) = /p,0X 



(4) 



where, denotes the number of atom x per unit volume, fp^Q is the planet occurrence 
rate at Solar metaUicity, and fp{[Fe/ H]) is the planet occurrence rate as a function of the 
metaUicity of the host star. We call this Functionl. To estimate the planet occurrence rate 
in NGC 6791 we also consider 



fpi[Fe/H]) 



fp,Q ^ 



Nh 

/p,0, for [Fe/H] < 



NFe \ 



, for [Fe/H] ^ 



(5) 



which we call Function2. Here we assume that the planet occurrence rate w ith respect to 



Fischer &: Valenti 



(120051) only for Solar 



the metaUicity increase as a power-law as found by 
or higher metallicities and is flat below Solar metallicities. For both cases the average planet 
occurrence rate is calculated considering all Kepler targets using /p^Kic = jiJ2fp{[Fe/H]). 
The expe cted planet occurrence rate /p,NGC 679i is also calculated for NGC 6791 using [Fe/H] 
= +0.30 ( iBoesgaard et al.ll2009l ). The boost of planet occurrence rate for giant planets larger 
than Saturn is then simply /p,NGC 679i//p,Kic- We find that this ratio is 4.9 and 3.1 assuming 
Functionl and Function 2, respectively. 

If the super-solar metaUicity boosts the planet occurrence rate for all planet sizes, then 
the total number of detectable transits will simply scale with the numbers presented in Table 
131 If the planet occurrence rate is boosted by the above factors in NGC 6791 only for giant 
planets (i.e., larger than Saturn), the expected number of transit detections will increase 
by a smaller factor depending on the number of giant planets, which in turn depends on 
the planet-size distribution. Table HI shows the expected numbers of transit detections by 
Kepler boosted due to the super-solar metaUicity of NGC 6791 assuming that the planet 
occurrence rate is boosted only for giant planets larger than Saturn and that the planet-size 
distribution is described by a power-law with exponent a = —1.99. 
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Table 4: The expected number of Kepler-detectable tran- 
siting planets if the planet occurrence rate in NGC 6791 
is boosted due to its high metallicity. We consider two 
different relations relating the giant planet occurrence 
rate and the metallicity of the host star (Equations H] 
andE]) denoted as Functionl and 2 respectively. For each 
relation we estimate the number of expected transit de- 
tections for various observation lengths assuming that 
the occurrence of planets more massive than Saturn is 
boosted due to metallicity. tobs denotes the duration of 
data collection by Kepler. For each of these data collec- 
tion durations the number of transit detections are shown 
both with and without the metallicity boost (assuming 
a = —1.99; Table [3]). All of the above numbers are given 
for the full cluster ("All Cluster"), the subregion of the 
cluster where pixel-level data is available (superaperture, 
a 200 X 200 pixels box centered on the cluster center; 
"Downloaded Cluster"; see Section [5731) . and the subre- 
gion of the cluster within the superaperture and outside 
'T' ^ ''"EE95, where the number of stars per EE95 is less 
than 1 ("r ^rEE95")- 



^obs (yr) 


Relation 


All Cluster 




Downloaded Cluster 
















r ^ rEE95 






Unboosted 


Boosted 


Unboosted 


Boosted 


Unboosted 


Boostec 


1 


Functionl 


10 ±2 


21 ±7 


8±1 


18 ±2 


3± 1 


6 ±2 


3.5 


(Equation H]) 


15 ±2 


31 ±7 


11 ± 1 


22 ±2 


4±2 


8 ±2 


8 




18 ±3 


34 ±7 


14 ±2 


26 ±3 


5±2 


9 ±2 


1 


Function2 


10 ±2 


16 ±5 


8±1 


14 ±2 


3± 1 


4±2 


3.5 


(Equation [5]) 


15 ±2 


23 ±5 


11 ± 1 


17 ±2 


4±2 


6±2 


8 




18 ±3 


26 ±6 


14 ±2 


20 ±3 


5±2 


6±3 
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Figure 10. Left Panel: Two dimensional image of our simulated model of NGC 6791. The grey dots denote all cluster objects 
(single and binary stars). The red circles denote planet hosts detectable {SNR > 7) by Kepler. Right panel: Magenta stars 
denote the positions of the individual stars for which the positions (within 23 arcseconds from the center of NGC 6791), 
magnitudes {Kp > 16), and the surface gravities (logg > 3.5) are consistent with their being main-sequence stars in NGC 6791, 
and for which data is being (or will be) downloaded by Kepler. For both panels the dashed square denotes the 200 X 200 pixel 
block of super apertures where pixel level data is being downloaded. The dashed circles from small to large are the rc, r^i, 
'"EEOSi and the tidal radius (partly visible in the scale) for reference. Note that the half-light radius is just inside the radius 
where the number of stars per EE95 is < 1. All lengths are denoted in units of Kepler pixels 3.98". 



5.3 Effects for finite data download by Kepler 

Our analysis in Section [S] shows that the best chance for Kepler to detect transits in 
NGC 6791 is around bright MS stars (16.5 < Kp < 20; Figure [7]). All predicted numbers of 
transit detections from our simulations in Sections 15.11 - [5l2] have been calculated assuming 
that data for all stars in NGC 6791 up to Kp = 20 is available. In practice, not all needed 
pixels are downloaded. Kepler has been observing a few regions in the FOV through large, 
custom apertures where pixel level data is downloaded. These are called "superapertures" 



(http://keplergo.arc.nasa.gov/Blog.shtml). In particular for NGC 6791 a supe raperture of 
200 X 200 pixels block centered on the cluster center is being downloaded (e.g.. 



Stello et al 



20111 ). In addition to this superaperture data, data from about 400 individual stars consis- 
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tent to be MS stars in NGC 6791 (i.e., within the tidal radius = 23.1' of the center of 
NGC 6791, magnitude Kp > 16, and surface gravity \ogg > 3.5) are being downloaded by 
Kepler. Most of these individual stars with downloaded data are already included in the 
superaperture (Figure [TU]). The right panel of Figure ITU] shows the superaperture block, and 
the individual stars consistent to be MS stars in NGC 6791 for which data is available (or 
will become available as part of a GO proposal). Important cluster structural radii [vc = 3.3', 
Thi = 4:A', rEE95 = 6' and rt = 23.1') are also shown. The left panel shows a two dimensional 
image of our simulated model of NGC 6791 over-plotted with the superaperture block, and 
the above mentioned cluster radii. Detectable transits by Kepler will be around bright main- 
sequence dwarf stars in NGC 6791. Thus at the given age of NGC 6791 these host stars are 
among the lower-mass sub-population of the cluster (M^ ^ 1.2 M0). Mass-segregation dis- 
tributes these relatively lower mass stars at larger cluster centric distances compared to the 
stellar binaries and more-massive (and evolved) single stars in NGC 6791. Thus while red 
giants are predominantly near the center of the cluster, a large fraction 60%) of transit 
detectable planet host stars are outside Vhi- Ideally, the best place to search for these planets 
is between rEE95 and r^. However, due to the limited data availability this cannot be done. 
We find that between teeqs = 4.8' - 6' and the superaperture box 6.7' along a side of the 
square and 9.4' along a diagonal) 3 ± 1 transits could be detectable by Kepler with one year 
of data collection assuming a = —1.99 and no metallicity boost. Table H] lists the expected 
numbers (with and without metallicity boost) of transit detections for different observations 
lengths in different cluster regions (including the region between rEE95 and the superaperture 
box). However, note that most of the stars in the superaperture are not analyzed as part of 
the regular planet search pipeline. 

6 DISCUSSION 

We revisit the possibility of detecting planets around normal MS stars in clusters in the 
light of the highly successful Kepler mission. In particular, we focus on NGC 6791, an old, 
massive, and metal-rich open cluster. The fact that NGC 6791 is already being observed by 
Kepler makes it an ideal candidate to test hypotheses developed to explain previous planet 
searches that report a dearth of planets around normal MS stars in clusters. 

Using ~ 200 detailed numerical simulations of a cluster's evolution and planet-harboring 
stars we study effects of stellar encounters in rich open clusters. We choose a model that best 
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Figure 11. Cumulative number of observed Kepler candidates Nai,i,{< r) vs the cluster-centric positions in units of tidal radii 
(rt). Top panels show results for the region around NGC 6791 only. The bottom panels show combined results including all 
four open clusters in the Kepler FOV. The left, middle, and right panels include all stars with Kp < 20 in the Kepler FOV, a 
subset of that with contamination values c < 0.4, and another subset that are between the single-star a nd equal-mass binar y 
MSs on a T^ff vs. Kp diagram for each cluster, respectively. Black points show Kepler planet candidates (Batalha c t al.ir2012h . 
Error bars are standard la Poisson errors. The histograms in each panel show the expected number of planets Kepler could 
detect based on our models. These histograms are calculated from our models using a = —1.99 (Section 15.11 1. no metallicity 
boost due to the super-solar metallicity of NGC 6791 (Section 15.211 , and assuming that data for all stars up to Kp = 20 in 
NGC 6791 is available (Section 15.311 . The thick lines include only detectable transiting planets from our simulated models of 
NGC 6791. Thin lines include planets both in our simulated models and those assumed to be around fore ground/background 
stars based on the stellar densities inferred from the Besangon model of the Milky Way l lRobin et al.|[200^ ) . The black (solid) , 
red (dashed), and blue (dotted) lines show results based on 1, 3.5, and Syr observations by Kepler. The green (dash-dot) lines 
show results if the cluster had no planets. The vertical dotted lines show where the tidal boundary for NGC 6791 is to aid 
the eye. The two dashed vertical lines in each of the top panels show the positions of the boundaries along a side and along a 
diagonal of the square 200 X 200 pixels superaperture box centered on the center of NGC 6791 (Section l5.3llStello et al.ll201lf) . 



matches the observed properties of NGC 6791 from this large set of simulations performed 
on a broad multidimensional grid of initial conditions (Section [3], Table [1], Figure [1]). We find 
that planetary orbits are rarely disrupted solely via stellar encounters in open clusters for a 
large range in cluster mass and stellar density (Table [2]). For clusters similar to NGC 6791, 
only about 10% of planetary orbits with relatively large semimajor axes (op ^ 10 AU) are 
likely to be excited via strong stellar encounters (FigureH]). A small fraction (~ 1%) of large- 
Op orbits may get ionized (Table [2]). However, the bulk of the planetary orbits, especially the 
small-Op orbits (op < 1 AU), detectable via transit searches, remain undisturbed (Figures [3] 
andH. 
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We find that the number of free-floating planets can range from a few to ~ 100 in a cluster 
like NGC 6791 depending on the age of the cluster. In our best-match model of NGC 6791 
the highest fraction of free-floating planets bound to the cluster potential, fpjf abound, grows 
from an initial value of zero by construction to the maximum value of 0.4 (Tabled]). However, 
depending on the cluster properties and age fpjf, bound can be higher (Figure [5l Table [2]). The 
total number of planets in a cluster steadily decreases with the age of the cluster mainly due 
to Galactic tidal stripping of the host stars (Figure |6]). These planet-host stars then leave the 
cluster to populate the field. Planetary orbits with a > 100 AU (and up to 1000 AU in our 
simulations) can remain bound to the host stars in a cluster similar to NGC 6791 (Figure 

ED- 



6.1 Prospects for cluster-planet detection from Kepler 

We estimate the expected number of planets Kepler may discover in NGC 6791 using the 
best-match initial conditions obtained from our grid of simulations to create models of 
NGC 6791 with planet frequency and planetary propertie s (Section l2.2p guided by the cur- 
rent Kepler observations in the field ( iBatalha et al. 1 120121 ). If the planet occurrence rate in 
NGC 6791 is the same as observed in the field by Kepler, then the existing Kepler data could 
detect transits of between a few to 10 planets (Figure [TJ Table E]) depending on the intrinsic 
distribution for planet sizes analyzing 1 year of data. An extended Kepler mission over 8 yr 
(or more) is projected to increase this yield by a factor of 2. However, due to the old age 
(8 ± 1 Gyr) and large distance (4kpc) of NGC 6791, most of these detectable-planet-host 
stars are relatively faint {Kp ^ 16.5, the MS turn-off for NGC 6791; Figure [2]). Hence, to 
attain such an yield, Kepler data for relatively faint stars {Kp < 20) must be analyzed prop- 
erly. These numbers are expected to increase due to the high metallicity of NGC 6791 (given 
the positive correla tion between the intrins ic planet occurrence rate and the metallicity of 



the host stars, e.g.. 



Fischer &: Valenti 



20051 ) (Table H]). On the other hand, since data is not 



available for all stars in NGC 6791 the actual number of transit detection will be reduced 
(Table m). 

Although a few giant planets may be detectable around low-luminosity giant stars, most 
planets should be detected around MS stars close to the MS turn-off (Figure [7j). When using 
only 1 yr observation by Kepler a large fraction (35%) of the detected planets will likely 
be gas giants. With longer observation times the expected number of detections for smaller 
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planets grows and the median Rp of detected planets decreases from 7 i?® after 1 yr to about 
4i?0 after Syr (for a planet-size distribution with a = —1.99; Table With an extended 
mission of 8 yr a few transits of planets as small as -Rp ~ 2 could be detected in NGC 6791 
(Figure [H]). 

About 70% of the detectable planets are expected to reside outside the core (rc ~ 3.3') 
of NGC 6791, and a little above 60% of the detectable planets are expected to reside outside 
the half-hght radius {r^i ^ 4.4') of NGC 6791 (Figure ED • Between 30-50% of all otherwise 
detectable planets reside outside the 2D projected distance teeqs = 4.8' - 6' from the cluster 
center beyond which the number of stars per EE95 is less than 1. Hence, crowding should 
not be a limiting factor in the detection of planets in NGC 6791 using Kepler, once the data 
between rEE95 and the superaperture boundaries (6.7' along a side and 9.4' along a diagonal) 
is fully analyzed. 

Note that the MS planet hosts, are among the lower mass (M^ ^ 1.2 Mq) subpopulation 
in the cluster given the old age of NGC 6791. Hence mass segregation distributes them at 
relatively larger separations from the cluster center compared to the stellar binaries and 
more massive evolved stars (Section |5]). Ideally, the best place to search for transits in this 
cluster would be between teeqs (to avoid crowding) and ~ O.Sr^ (since there are simply not 
many stars beyond this radius). Hence, based on our results we suggest the following. 

1. Relatively fainter stars {Kp < 20) should also be properly analyzed. 

2. Ongoing observations of the clusters in the Kepler FOV, especially the relatively massive 
clusters NGC 6791 and NGC 6819 are maintained not only for astrophysics, but also to 
search for exoplanets. 

3. The pixel- level data obtained from the superapertures are analyzed for existence of tran- 
siting planets. Especially, data for regions outside teeqs where there is reduced crowding. 

4. For detection of exoplanets it is more important to get data from the outer (e.g., outside 
the Thi) regions of a cluster. Considering the limitation of the number of pixels that can be 
downloaded, we suggest that some pixels from the center of the superapertures be reassigned 
to search for planets around individual MS stars likely to be cluster members and located 
outside the superaperture block. 

Kepler results can determine whether planets in short period orbits are as common 
around cluster stars as they are around field stars. Even a null result from Kepler would 
be interesting since that would potentially provide stronger constraints than the existing 
results on the frequency of planet occurrence around MS stars in star clusters. 
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6.2 How likely is a planet discovered near a cluster an actual cluster member? 



In Figure [TT] we show a preliminary analysis of the Kepler planet candidates near star 
clusters as a function of their sky positions relative to a cluster. Planet candidates are from 
(IBatalha et al.ll2012l ) and are based on searching Kepler data from Q1-Q6. The top panels 
show the region around NGC 6791 only, and the bottom panels combine the regions around 
four open clusters (NGC 6866, NGC 6811, NGC 6819, and NGC 6791) in the Kepler FOV. 
We plot the number of observed Kepler candidates (A^obs) as a function of the projected 
distance from the cluster center (s). The cluster-centric distances are given in units of the 
cluster's tidal radius(ii) (rj). Cluster parameters, inc luding positions, tidal radii, distances. 



ages, and redennings are taken from the literature (jKalira 



2005 



Hole et al 



2009 



Piskunov et al. 



2008 



Platais et al 



et al. 



2011 



2001 



Kharchenko et al. 



WEBDA) The observed 



points are drawn for three sets of target stars defined in the following way. In the left panels 
we include all stars with Kp < 20 that were observed by Kepler. In the middle panels we 
limit the sample to include only stars with contamination c < 0.4 to reduce the effect that 
stray light from nearby stars may have on the frequency of Kepler candidates, which may be 
significant in clusters like NGC 6791. This value of 0.4 is chosen fairly arbitrarily, and a more 
thorough analysis is needed to determine the proper cutoff in c to draw robust conclusions. 
Finally, in the right panels, we limit the Kepler sample to only include stars that would lie 
between the s i ngle-s tar and equal-mass binary sequences of the respective clusters. We use 



Marigo et al 



(|2008[ ) isochrones in log(Te//) vs. Kp space to perform this selection. 



Although the vast majority of the Kepler observed stars (and especially those targeted 
for the "EX" program) have Kp < 16, we include all stars with Kp < 20 observed by 
Kepler to remain consistent with our theoretical expectation that planets have the highest 
detection probability between 16 < Kp < 20 in NGC 6791. Given the old age and distance of 
NGC 6791, the cluster MS begins at Kp ^ 16 (Figure |2]). Consequently, the majority of the 
observed Kepler candidates in Figure [TT] for NGC 6791 must be giants, or possibly subgiants, 
if they are really around cluster members. Our numerical models suggest that it is unlikely 
for Kepler to detect more than about 3 planets around giant host stars in NGC 6791 (Figure 
[7j). Kepler has not yet reported strong planet candidates transiting stars that are within Irt 
from the center of NGC 6791 and have magnitudes and temperatures consistent with a MS 
cluster member (top- right panel of Figure [TT]) . Given that our models predict significantly 
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more planet detections possible by Kepler for the MS stars in NGC 6791, we suggest that a 
detailed analysis of these fainter stars for planet transits would be highly valuable. 

If a transit is detected in the direction of a cluster, it is not a priori clear whether the 
planet is indeed around a cluster star or around a foreground/background star. We find that 
if a transiting planet is detected within 0.5 rt from the center of NGC 6791 and the host star 
has properties such that it will lie within the single and equal-mass binary MSs in the CMD 
of NGC 6791, then the planet is very likely to be around a true cluster member. About 60% 
of all transiting planet detections around stars satisfying this criteria and within rt from the 
center of NGC 6791 are expected to be around true cluster members (Figure HH top-right). 
Given that there are about 10 Kepler candidates within the tidal radii of the 4 open clusters 
in the Kepler FOV (Figure [TTl bottom-right) around stars with properties satisfying the 
CMD based limits described above, about 6 of these systems may be true cluster members. 
The relative contributions from the cluster stars and the foreground/background stars do 
not change significantly if the power-law exponent for the planet-size distribution is changed 
in the range explored in this study (Table |3]). However, note that this part of the analysis 
has a few caveats as described below. 

We compute the predicted number of Kepler detectable planets transiting foreground/background 
stars superimposed on our model of NGC 6791 (thin lines in the top panels of Figure [TTl) . 
To produce these curves, we embed the simulated cluster within a uniform field of back- 
ground/foreground stars and treat these data in the sar ne manner as for t he Kepler obser- 



vations. We use the Besangon model of the Milky Way ( iRobin et al 



20031) to estimate the 



i^p-dependent surface densities of s t ars in the direction of NGC 6791. Using the transfor- 



mation equations from 



Smith et al 



( 120021 ) we convert the Johnson B, V magnitudes given 
in the Besangon model to SDSS g, r, and then use the standard Kepler conversions to find 
Kp. We estimate the surface densities of background/foreground stars in the direction of 
NGC 6791 as a function of Kp. The surface density of the expected Kepler detections for 
planets transiting these background/foreground stars per solid angle is then estimated using 
the same method described in Section |2] for Kp values between 5 < Kp < 20. These esti- 
mated surface densities are then summed to obtain the integrated surface density of planet 
detections for non-cluster stars over the full range in 5 < Kp < 20. The expected cumulative 
number of transit detections for non cluster stars as a function of the distance (r) from the 
cluster center is then found simply by multiplying the expected number of transit detections 
per sky-projected area and vrr^. 
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We find that the total number of background/foreground stars estimated from the Be- 
sangon model can be lower by ~ 30% than that counted from the Kepler Input Catalog 
(KIC) by centering at the same Galactic latitude but away from the cluster. Even the num- 
bers of stars estimated from the KIC catalog using a center at different nearby locations 
at the same Galactic latitude can vary by about 10%. These differences in the number of 
stars are directly translated to differences in the estimated number of detectable transits 
around non cluster members. For example, using the KIC stars within an angular radius 
of 23. l' centered around a point away from NGC 6791 but at the same Galactic latitude 
as NGC 6791, this model would estimate an expected number of detectable planets within 
r < 1 r^, transiting non cluster stars with Kp < 17 to be 49. Using the same cut-offs for Kp 
and r, but using the Besangon model the estimated number of detectable transiting planets 
is 34. We caution that, both of the above estimates are likely upper limits on the number of 
detectable planets transiting non-cluster stars, since it assumes that all non-cluster stars are 
dwarfs. A significant fraction of the foreground/background stars will be giant or subgiant 
stars for which the minimum detectable planet size is significantly larger than for a MS star 
of the same magnitude. Thus, the thin lines should be treated as rough estimates for an 
upper limit to the number of contaminating planets around foreground/background stars. 
Fortunately, even these results suggest that planet candidates within 0.5 are likely cluster 
members if the selected sample has properties consistent with their being on the MS of the 
cluster. 

A more complete and careful study of the frequency of planets in the fields near the open 
clusters in the Kepler survey is of great interest (but beyond the scope of this paper). We 
stress here that a complete study of any putative trend between planet occurrence rate and 
position relative to a cluster must be mindful of a number of caveats (in addition to those 
mentioned above) related to how the sample of target stars are chosen before drawing secure 
conclusions. For example, the search for planets in clusters is still incomplete (for NGC 6791 
it has barely started) due to technical details related to targeting and the pipeline. A careful 
study must include proper cuts in the Kp values to fairly compare cluster stars with those 
in the field, and take into account differences in magnitude, stellar radii, contamination, 
pipeline selection effects between cluster and field stars, and preferably kinematic cluster 
membership information. 
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